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An ensemble with random n-body interactions is investigated in the presence of symmetries. A 
striking emergence of regularities in spectra, ground state spins and isospins is discovered in both 
odd and even-particle systems. Various types of correlations from pairing to spectral sequences and 
correlations across different masses are explored. A search for interpretation is presented. 
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Recent progress in ab-initio treatment of light nuclei 
show unambiguously the essential role played by the 
three and four-body forces [3, [2] ■ With the advancement 
of the mesoscopic physics and with the ability to manip- 
ulate interactions the question how n-body interactions 
play out in the many-body physics becomes increasingly 
important. 



It is known since the Wigner-Dyson random matrix 
theory (RMT) 0| , see also [j, 0| , that complex configu- 
ration mixing driven by many-body forces have generic 
statistical features which also depend on the nature of 
interaction Symmetries have a robust effect, which is 
not fully understood. Recently, a remarkable geometrical 
ordering have been numerically found in the Two-Body 
Random Ensembles (2-BRE)[7j]. The main question is 
centered around the disproportionally large probability 
for the ground state (g.s.) spin J to be zero, Jo = 0, in 
the even-particle system. Manifestations of the symme- 
try were seen in features of the mean-field @, Q| , and in 
transition probabilities indicating vibrational and rota- 
tional low- lying spectra [T3|. Starting from the pioneer- 
ing paper [7| over a hundred works were published by 
different groups striving to understand the emergence of 
symmetries, their role and origin in the 2-BRE. The sum- 
mary of these efforts may be found in reviews [13, 0, E2] . 
The success is mixed, we understand that pairing [13j and 
time-reversal are not the primary causes. The boson 
approximation of fermion pairs and chaos arising from 
complex geometrical couplings, the geometric chaotic- 
ity, provide some qualitative understanding of the trend; 
more quantitative explanations have been suggested in 
Refs. [H, [lB| via numerical observations of geometrical 
correlations and diagonalization of the individual inter- 
action terms. Nevertheless, the simple question of sym- 
metry and chaos asks for a simple answer which is still 
missing. Driven by the quest for understanding of in- 
terplay between symmetry and many-body complexity 
and the general interest and importance of many-body 
interactions in modern physics in this work we address 
the n-body Random Ensembles (n— BRE) with symme- 
tries. The evolution of spectrum and level spacing for 
the Gaussian n-BRE without symmetries is discussed in 
Ref. 0. 
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where operators T^m ( a ) are n-body creation operators 
coupled to a total angular momentum L and magnetic 
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projection M, (a) 
here 1 is the single-particle index. For n = 2 the coeffi- 
cients Cf 2 M are proportional to the Clebsch-Gordan coef- 
ficients and index a is uniquely defined by single-particle 
levels involved. For n > 2 the index a must include ad- 
ditional information about the coupling scheme. Here 
we define the basis set of normalized n-body operators 
Tlm ( a ) fr° m a full se t of orthogonal n-body eigenstates 
Tlm (oO|0) of an arbitrarily chosen reference two-body 

(2) 

Hamiltonian Hq , solved numerically. The exact form 
of this Hamiltonian is irrelevant as long as it preserves 
rotational and other symmetries of the problem. 

The n-BRE is defined as a set of n-body interaction 
Hamiltonians that for an n-particle system leads to a 
Gaussian Orthogonal Ensemble (GOE) within every sym- 
metry class. The interaction strengths V£ (a, 0) are 
selected at random with normal distribution centered 



at zero 
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and unit diagonal variance 
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time reversal symmetry sets 
thus we assume ordered labels a 
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is Hq independent. The variance defines the energy 
unit. Our ensemble extends the n-body Embedded GOE 
[l8| by incorporating symmetries. Inclusion of the 
isospin symmetry is straightforward. 

Through the text we use P [event] to denote chances of 
observing a certain event in the n-BRE. Most commonly 
we discuss J(N)q which is an event where g.s. of an 
iV-particle system has spin J; the subscript reflects the 
order in excitation energy, 0-for g.s. 1-first excited state, 
and etc. Where obvious we omit the explicit reference to 
the particle number N. 

We start with a single j-level model of 2j + 1 degen- 
eracy with N identical fermions. In Fig.Q] some exam- 
ples are shown. In the special case of n = 1, the mean- 
field only, all many-body states are degenerate. Here 
by definition we assume an equiprobable ordering of the 
states in the spectrum. Thus, the probability to observe 
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a g.s. with a spin J, P [Jo] is by definition proportional 
to the number of many-body states with this spin in the 
model space d(J). Following the semi-classical consider- 
ation of the random- walk- type vector couplings, see also 
[Hi d(J) ~ (2J + 1) exp[-3 J( J + l)/2Nj(j + 1)] which 
disfavors both J = and the maximum possible spin 
Jmax = N(2j + 1 — N)/2. The equiprobable ordering is 
in drastic contrast to the 2-BRE where g.s. is most likely 
to be Jo = and the probability for Jo = J max is large, 
Fig.ffl The preponderance of P [0 ] becomes stronger 
for the n-BRE with higher n, and at the same time the 
chances of J max as g.s. quickly diminish. For n — 4 and 
5, apart from the dominant Jo = 0, the states with spins 
Jo = 2,4,6,8 have small, few percent-level, chances to 
appear as g.s. 

The middle panel in Fig.Q] corresponds to the TV = 7 
odd-particle system. The lowest bar in the stacked his- 
togram shows the P\jo] (where j = 19/2) that can be 
interpreted as a single particle coupled to the J = 
TV = 6 core. While for the 2-BRE P [19/2 ] = 12.4(4)% 
(note that P [13/2 ] =14.5(4)%), the P [19/2 ] is big- 
ger for n = 3, 4 and 5. The maximum spin probability 
F[(91/2) ] is enhanced for the 2-BRE but similarly to 
the even-particle system declines for larger n. Here, and 
on some occasions below we include statistical errors due 
to a limited number of random realizations tested. 

An important case of n = TV is not show in Fig.Q] be- 
cause of its drastic contrast to the n < TV situations. 
Here, within each symmetry class the Hamiltonian ma- 
trix is represented by the GOE for which with increas- 
ing dimensionality d{J) the distribution of eigenvalues 
quickly conve rges to the Wigner semicircle with the ra- 
dius y/2d(J). A detailed quantitative analysis can be 
done using the RMT but it is clear that the probability 
P [Jo] strongly favors those spins J with the highest di- 
mensionalities d(J). For example, for TV — 8 j — 19/2 the 
highest dimensions are 179, 178, 173, and 169 for spins 
J =12, 14, 16, and 10; the corresponding probabilities 
P [J ] in the 8-BRE are 41, 35, 13 and 5%, respectively. 
For lower n = TV — 1 the change in the spin statistics of 
g.s. occurs abruptly. 

Although pairing was ruled out as an explanation for 
the 2-BRE given its dominance in realistic systems it is 
worth revisiting this question for the n-BRE. In Fig. [2] 
the evolution of seniority s (the number of unpaired par- 
ticles) as a function of n is shown. The seniority s is 
evaluated using expectation value in the state of interest 
|TV, a) of the L — pair operator, which is unique in a 

single-j (N,a\T$ 1 T$\N,a) = (TV — s)(2j + 3 - TV - 
s)/{2(2j + l)}. 

Surprisingly, the paired state is favored by the three- 
body forces. Significant lowering of the seniority for 
n = 3 is a robust result in all cases considered Fig. [21 
Given that the P[0q] goes up for higher n while s is in- 
creasing for n > 3 pairing still does not explain the pre- 
ponderance of the zero spin g.s. The structure of the two- 
body forces on a single-j level is known to facilitate se- 
niority conservation, only about a third of the j — 1/2 in- 
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Figure 1: (Color online) The stacked-bar histogram showing 
the distribution of ground states with a given angular mo- 
mentum J for a single j-system, with j and TV as marked. 
The histograms are shown as a function of n for different n- 
BRE. For the n = 1 case we assign probability proportionally 
to the number of states with that spin. The stacking order 
is marked, and corresponds to the increasing J starting from 
Jmin , with the exception of the odd- TV where Jmin = j- The 
spins J with P [Jo] < 2% are not separately identified, their 
cumulative probability is shown with the white bar, labeled 
"other". 
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Figure 2: The ensemble averaged seniority s is shown as 
a function of n for the n-BRE. Left panel corresponds to 
j — 19/2 with connected by line points for different parti- 
cle numbers. Only (J m in) ground state realizations are se- 
lected, namely Jo = for N even and Jo = 19/2 for odd. On 
the right the half-occupied system TV = 8 j = 15/2 is exam- 
ined, which shows the average seniority in realizations with 
g.s Jo = 0, 2 and 4. For the n = 1 case, (equiprobable ground 
state distribution) the average seniority of states with a given 
spin is quoted. 



dependent linear combinations of interaction parameters 

(2) 

V L mix seniorities, seniority is a g ood quantum number 
for any interaction on j < 9/2 [20(|. There is lowering of 
seniority against the average (n = 1) for n = 2 and per- 
haps explanation for n = 3 pairing enhancement lies in a 
similar seniority conserving structure of the three-body 
interactions. This question is a subject for future work. 

The particle number TV and the Casimir operators for 
the symmetry groups such as J 2 and T 2 are conserved 
quantities, and in each random realization of the Hamil- 
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tonian are expected to describe its coherent part. In the 
2-BRE the particle-hole transformation can be used to 
obtain a coherent J 2 -dependent component of the Hamil- 
tonian {H^ 2 >)j = V\ J 2 . The moment of inertia coefficient 
Vi is given via angular momentum recoupling coefficients, 
for a single-j it is 



3(2£+l) 
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The 2-BRE g.s. spin systematics follows from here; see 
also geometric chaoticity arguments The Vi, be- 

ing a sum of random numbers with normal distribution 
centered at zero, is itself a normally distributed random 
variable with equal chances of being positive and nega- 
tive. Thus the P[0 ] = P[(J ma x)o] = 1/2. This predic- 
tion qualitatively describes the observations, although it 
is distorted by other incoherent interaction terms . 

The emergence of the moment of inertia in each real- 
ization of interaction is further supported by the mass 
number independence of the 2-BRE results. Indeed, the 
moment of inertia ^ is particle-number independent, 
thus for a given realization of the two-body Hamilto- 
nian if V\ > the g.s. spin is zero J(N)q = for any 
even TV (j for odd); similarly for V% < the g.s. is 
J max (TV) . To quantify the correlation between the statis- 
tics for different TV we consider the joint probability that 
all systems from TV — 5 to half-occupied, simultane- 
ously have the maximum spin P [J ma x (5)o, J ma x(6)o, ■■■}■ 
The same can be done for the minimum spin J m in 
which we define to be zero for an even TV and s.p. j 
for odd. Because of the particle-hole symmetry the 
eigenvalues of the 2j + 1 — TV system apart from a 
monopole constant shift in energy are identical to the 
TV-particle system. In an uncorrelated case the joint 
probability is a product of the independent probabilities 
P [J(/Vi) , J(TV 2 ) , . . . } = n« P \J( N i)o] which is gener- 
ally small. On the other hand, if strong correlations 
exist the joint probability is of the order of individual 
P[J Q }. In our example j = 19/2 with TV = 5, 6... 10 
P[J m ax(5) .--J,nax{lO) } = 6.6% in the 2-BRE, while 
uncorrelated product is four orders of magnitude smaller 

P [J m ax(5) ] P [Jm Q *(6)o] ...P [Jm«x(10)o] = 2.1 -ltT^. 

In Tab. [J we show the total weighted correlation which for 
a general set of events J, J', J", J"' ... is defined as 



C[J,J',J",...] 



log(P[J]P[J']P[J"]...) 

log (P[J,J',J",...}) 



1. (2) 



If ground states in k systems with different masses are 
not correlated then C — 0. While in the case of full 
correlation the ratio of logarithms is approximately equal 
to k and C = k — 1, meaning that the g.s spin in one 
system is sufficient to predict g.s. spins in all other k — 1 
different mass-systems. 

The situation with n-body forces is more complex as 
coherent higher order terms appear. For the 3-BRE 
(HW)j = (V'N + V)J 2 and the moment inertia is 
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16.0 


10.1 


36.3 


2.9 


7.7 


0.2 


6 


52.3 


10.5 


66.4 


3.1 


83.0 


0.0 


7 


12.4 


11.8 


39.1 


4.8 


33.0 


0.5 


8 


42.7 


12.1 


63.2 


5.0 


84.3 


1.1 


9 


9.5 


12.3 


31.1 


6.5 


33.7 


2.3 


10 


31.2 


11.5 


48.6 


7.1 


65.5 


2.6 


C 


1.883 


3.806 


1.560 


3.266 


0.435 


0.000 



Table I: Summary of g.s. statistics for minimum and maxi- 
mum spin, and correlations across different mass numbers TV 
for a single j = 19/2 valence space with 2,3, and 4-body ran- 
dom interactions. The first six rows correspond to P[J(N)] 
expressed in percent. The lowest row shows correlation C, for 
all particle numbers TV = 5, 6, 7, 8, 9, and 10. Columns reflect 
the type of ensemble and J m in or Jmax ■ 



particle-number dependent. However, the 3-BRE is ex- 
pected to be similar to the 2-BRE where the g.s. spin 
statistics is dictated primarily by the sign of the mo- 
ment of inertia which equally favors both J m in and J ma x ■ 
Although the effect is reduced by incoherent interac- 
tion components the preponderance of J m in and Jmax 
is seen m Fig.H and Tab.[U The TV-depend ence reduces 
the amount of correlation between systems of different 
masses, see Tab. HI However, the variations of TV which 
is relatively large and positive are usually insufficient to 
change the sign in the moment of inertia and therefore 
correlations between systems of different TV are strong. 

The angular momentum dependent part in the 4-body 
forces is given by a more complicated expression with 
four possible, normally distributed, interaction terms 
(HW)j = (V{'N 2 +V{N + Vi)J 2 + V 2 J 4 - The coeffi- 

cients V2,Vi,V{, and V\ are given by correlated nor- 
mal distributions but because of geometric complexity 
in the recoupling coefficients and a large number of 4- 
body interaction parameters these correlations are small. 
The P[(Jmax)o\ is expected to be significantly reduced, 
with more chances going to the P [(J m in) ] which is al- 
ways a local minimum for non-negative J 2 , Tab.|l] and 
Fig.ffl Only for a large particle number the term V/'TV 2 J 2 
in the Hamiltonian seems to dominate enhancing the 
P [(Jmax) ] , Tab.|H Because of the strong TV-dependence 
in the moment of inertia and presence of the centrifugal 
distortion V 2 for n = 4 there is much less correlation be- 
tween g.s. spins for systems with different masses, Tab.|H 
The prevailing spin sequence of low-lying states pro- 
vides yet another information about the coherent sym- 
metry structure of interactions. The chances to find g.s. 
with spin Jo = followed by the first excited state with 
J\ = 2 and then by the second excited state J 2 — 4 are 
high. In Tab.Ull this probability P[0o,2i,42] is shown 
along with the average ratio of the excitation energies for 
the 2x and 4 2 states. The typical numbers between 2 and 
10 % reflect extremely high probability compared to the 
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E1/E2 
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E1/E2 


P 


E1/E2 


6 


3.7(2) 


0.55 


4.2(2)* 


0.69 


4.4(7)* 


0.69 


8 


4.2(2)* 


0.59 


5.5(2) 


0.67 


7.4(11) 


0.75 


10 


2.1(1)* 


0.72 


5.2(2) 


0.69 


7.3(10) 


0.62 



Table II: Probability of finding the three lowest states as a 
sequence 0,2,4, P[0o,2i,42] , labeled as P, expressed in per- 
cent, and the ratio of excitation energies between 2i and 42 
states. In all cases the sequence 0,2,4 in the most likely g.s. 
sequence except for those marked with *. 



chances of finding this sequence in a random list of spins. 
However, the term "sequence" must be used with reserva- 
tions, the weighted correlational entropy ^ between the 
joint P [0o, 2i, 4 2 ] and the product of independent P [0o], 
P[2i] and P [4 2 ] is only about C[0 ,2i,4 2 ] » 0.4 in all 
cases. The number is comparable to unity showing def- 
inite correlations, but they are not strong and the high 
probability to observe the sequence comes partially from 
the independently high chances of finding low-lying states 
with angular momenta 0, 2 and 4. These states do not 
always form a rotational band which would lead to the 
ratio of excitation energies E1/E2 = 0.3. It is likely that 
the members of the ground state band are often higher 
in energy and further work is needed to identify them. 

We conclude this work by showing the statistics of g.s. 
quantum numbers for systems with the rotational and 
isospin symmetries in Fig.[3l The preponderance of the 
symmetric g.s. with (JT) = (0 0) is robust for every 
particle number divisible by 4. Although geometrically 
(0 0) state is possible in j = 19/2 N — 10 it does not 
appear as g.s. The possibility of a-type correlations is 
to be investigated in the future. The g.s. in the N = 9 
system can be explained as a single particle coupled to 
the N = 8 core, thus Tq = 1/2 is favored along with 
(9/2 9/2) and (49/2 1/2), the maximum isospin or spin 
states. For N = 10 the P[(JT) ] is only substantial 
when either J = or T = with clear preference to 
the quantum numbers of two coupled identical j = 9/2 
fermions ( J + T is odd). The preponderance of the min- 
imum or maximum in either spin or isospin appears in 
all cases. Probability to find g.s. with J m ax or T max di- 
minish with higher n. The overall picture is in qualitative 
agreement with the hypothesis discussed above: that the 
powers of J 2 and T 2 operators appear with random sign 
as property of the coherent interaction components in 
each realization. The mixed terms, such as J 2 T 2 do not 
single out a particular ( J T) pair. In relation to a more 
complex shell model studies it is interesting to mention 
a p-shell (j — 3/2 and j — 1/2) odd-odd 10 B case, where 



results are similar to the above study favoring g.s. of cou- 
pled two-particle quantum numbers. For a degenerate p 
levels P [(1+0)q] = 19.1(4)%, P [(3+0) ] = 36.6(6)%, and 
P [(0+1) ] = 21.3(5)% for 2-BRE; the same numbers are 
25.6(5), 30.5(6), 24.5(5) for 3-BRE; and 30.3(6), 29.9(6), 
23.2(5) for 4-BRE, respectively. The chances of (J 7T T) 

j=9/2 N=10 
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Figure 3: Statistics of g.s. quantum numbers labeled as (JT) 
for isospin 1/2 fermions on a single j = 9/2 level. The panels 
from left to right correspond to N = 8, 9, and 10. The n = 1 
is equiprobable distribution. Notations are similar to Fig.[TJ 



being (l + 0) grow with n, although the statistics has little 
to do with the realization chosen by nature [l|, 0] • Fur- 
ther exploration of more complex models will be reported 
elsewhere. 

To summarize, in this work the ensembles of fermions 
interacting randomly with symmetry conserving many- 
body forces have been considered. The preponderance 
of symmetry dominated ground state was observed. The 
effect is generally stronger for the higher n-body forces. 
The possibility of correlated paired ground state is dis- 
cussed and a surprising enhancement of pairing with 
three-body forces is observed. An explanation based 
on the coherent components of interaction is suggested, 
which qualitatively describes the results. The strong cor- 
relations between systems of different particle-number, 
mass dependence of probabilities, and correlated se- 
quences of states support this theoretical hypothesis. 
This work opens many future avenues for investigation, 
studies and ideas @, 0, EE Ell from random two-body 
ensembles can be extended. The onset of coherence: pair- 
ing, in isovector or isoscalar form, a-particle four-body 
clustering, shape properties and vibrations are all inter- 
esting and important questions for future investigations. 

The author is grateful to V. Zelevinsky, D. Mulhall, 
and V. Abramkina for motivating discussions. Support 
from the U. S. Department of Energy, grant DE-FG02- 
92ER40750 is acknowledged. 



[1] P. Navratil and W. E. Ormand, Phys. Rev. Lett. 88, [2] S. C. Pieper, V. R. Pandharipande, R. B. Wiringa, and 
152502 (2002). J. Carlson, Phys. Rev. C 64, 014001 (2001). 



5 



[3] C. E. Porter, Statistical theories of spectra: fluctuations; 
a collection of reprints and original papers, with an in- 
troductory review (Academic Press, New York, 1965). 

[4] T. Guhr, A. Miiller-Groeling, and H. Weidenmiiller, 
Phys. Rep. 299, 189 (1998). 

[5] T. A. Brody, J. Flores, J. B. French, P. A. Mello, 
A. Pandey, and S. S. M. Wong, Rev. Mod. Phys. 53, 
385 (1981). 

[6] J. B. French and S. S. M. Wong, Phys. Lett. B, 35, 5 
(1971). 

[7] C. W. Johnson, G. F. Bertsch, and D. J. Dean, Phys. 

Rev. Lett. 80, 2749 (1998). 
[8] V. G. Zelevinsky, Nucl. Phys. A 555, 109 (1993). 
[9] R. Bijker and A. Frank, Phys. Rev. C 65, 044316 (2002). 
[10] V. Zelevinsky and A. Volya, Phys. Rep. 391, 311 (2004). 
[11] T. Papenbrock and H. A. Weidenmuller, Rev. Mod. Phys 
79, 997 (2007). 



[12] Y. M. Zhao, A. Arima, and N. Yoshinaga, Phys. Rep. 
400, 1 (2004). 

[13] D. Mulhall, A. Volya, and V. Zelevinsky, Phys. Rev. Lett. 

85, 4016 (2000). 
[14] R. Bijker, A. Frank, and S. Pittel, Phys. Rev. C 60, 

021302(R) (1999). 
[15] Y. M. Zhao, A. Arima, and N. Yoshinaga, Phys. Rev. C 

66, 034302 (2002). 
[16] T. Papenbrock and H. A. Weidenmuller, Phys. Rev. Lett. 

93, 132503 (2004). 
[17] K. K. Mon and J. B. French, Ann. Phys. 95, 90 (1975). 
[18] V. K. B. Kota, Phys.Rep. 347, 223 (2001). 
[19] H. A. Bethe, Phys.Rev. 50, 332 (1936). 
[20] A. Volya, Phys. Rev. C 65, 044311 (2002). 
[21] A. P. Zuker, Phys.Rev.Lett. 90, 042502 (2003). 



